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ABSTRACT 

We present a systematic search for parallax microlensing events among a total of 512 
microlensing candidates in the OGLE II database for the 1997-1999 seasons. We fit 
each microlensing candidate with both the standard microlensing model and also a 
parallax model that accounts for the Earth's motion around the Sun. We then search 
for the parallax signature by comparing the of the standard and parallax models. 
For the events which show a significant improvement, we further use the 'duration' of 
the event and the signal-to-noise ratio as criteria to separate true parallax events from 
other noisy microlensing events. We have discovered one convincing new candidate, 
sc33_4505, and seven other marginal cases. The convincing candidate (sc33_4505) is 
caused by a slow-moving, and likely low-mass, object, similar to other known parallax 
events. We found that irregular sampling and gaps between observing seasons hamper 
the recovery of parallax events. We have also searched for long-duration events that do 
not show parallax signatures. The lack of parallax effects in a microlensing event puts 
a lower-limit on the Einstein radius projected onto the observer plane, which in turn 
imposes a lower limit on the lens mass divided by the relative lens-source parallax. 
Most of the constraints are however quite weak. 

Key words: gravitational microlensing - galactic center 



1 INTRODUCTION 

Gravitational microlensing was originally proposed as a 
method of detecting compact dark matter objects in the 
Galactic halo (Paczyriski 1986). However, it also turned 
out to be an extremely useful method to study Galactic 
structure, mass functions of stars and potentially extra- 
solar planetary systems (for a review, see Paczyriski 1996). 
Most microlensing events are well described by the standard 
light curve (e.g. Paczyriski 1986). Unfortunately, from these 
light curves, one can only derive a single physical constraint, 
namely the Einstein radius crossing time, which involves the 
lens mass, various distance measures and relative velocity 
(see §3). This degeneracy means that the lens properties 
cannot be uniquely inferred, thus making the interpretation 
of the microlensing results ambiguous; this ambiguity is par- 
ticularly serious for the interpretation of microlensing events 
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toward the Large Magellanic Clouds (e.g. Sahu 1994; Zhao 
& Evans 2000; Alcock et al. 2000b, and references therein). 

Fortunately, some microlensing events are exotic, in 
that they deviate from the standard shape appropriate for 
a point source lensed by a single star (e.g. Paczyhski 1986). 
The deviations include binary microlensing events (e.g., Mao 
& Paczyhski 1991; Udalski et al. 2000; Alcock et al. 2000a), 
the finite source size events (Gould 1994; Nemiroff & Wick- 
ramasinghe 1994; Witt & Mao 1994) and parallax microlens- 
ing events. The importance of these non-standard events 
is that they allow one to derive extra constraints on the 
lens parameters. Parallax signatures in microlensing events 
arise when the event lasts long enough that the Earth's mo- 
tion can no longer be approximated as rectilinear during the 
event (Gould 1992; see also Refsdal 1966). Unlike the light 
curves for the standard events which are symmetric, these 
parallax events exhibit asymmetries in their light curves due 
to this motion of the Earth around the Sun. These events 
allow one to derive the physical dimension of the Einstein 
radius projected onto the observer plane (i.e., the solar sys- 
tem). The first parallax microlensing event was reported 
by the MACHO collaboration toward the Galactic bulge 
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(Alcock et al. 1995), and the second case (toward Carina) 
was discovered by the OGLE collaboration and reported 
in Mao (1999). Additional parallax microlensing candidates 
have been presented in a conference abstract (Bennett et al. 
1997). A spectacular microlensing event that exhibits par- 
allax signatures over a span of two years was reported by 
Soszynski et al. (2001). The MOA collaboration also dis- 
covered a parallax microlensing event toward the Galactic 
bulge (Bond et al. 2001). However, despite the importance 
of parallax events, there have been no reported systematic 
searches in the existing databases. This paper is an attempt 
to remedy this situation. We search for parallax signatures 
among 512 microlensing candidates discovered by Wozniak 
et al. (2001) using the difference image analysis. The outline 
of the paper is as follows. In §2 we briefly describe the mi- 
crolensing database, in §3 we describe our fitting and search 
procedures, and in §4 we describe our candidate parallax 
events, including both convincing and marginal cases. We 
also search the long-duration events that show no obvious 
parallax signatures and study their constraints on the lens 
parameters. Finally, in §5 we discuss the implications of our 
results and highlight observational issues in identifying par- 
allax events. 



2 DATABASE OF MICROLENSING EVENTS 

All observations presented in this paper were carried out 
during the second phase of the OGLE experiment with the 
1.3 m Warsaw telescope at the Las Campanas Observatory, 
Chile. The observatory is operated by the Carnegie Insti- 
tution of Washington. The telescope was equipped with 
the "first generation" camera with a SITe 3 2048x2048 
pixel CCD detector working in the drift-scan mode. The 
pixel size was 24^m, giving the scale of 0.417"pixel~'^. Ob- 
servations of the Galactic bulge fields were performed in 
the "medium" speed reading mode with the gain 7.1 c^ 
ADU~^ and readout noise about 6.3 e~. Details of the in- 
strumentation setup can be found in Udalski, Kubiak & 
Szymariski (1997). The majority of the OGLE-II frames 
were taken in the /-band, roughly 200-300 frames per field 
during observing seasons 1997-1999. A typical observing 
season for the Galactic bulge lasts between mid February 
to the end of October, which unfortunately produces gaps 
of more than a quarter of a year in length (for more de- 
tails, see Section 5). Udalski et al. (2000) gives full de- 
tails of the standard OGLE observing techniques, and the 
DoPhot photometry is available from the OGLE web site at 
http:/ / www, astrouw, edu.pl/ ~ ogle/ogle2/ ews/ ews. html. 

Wozniak et al. (2001) presented a sample of microlens- 
ing events from difi'ercnce image analysis of the three yeax 
OGLE-II bulge data. His difference image analysis pipeline 
is designed and tuned for the OGLE bulge data (Wozniak 
2000), and is based on the algorithm from Alard & Lup- 
ton (1998) and Alard (2000). The difference image anal- 
ysis pipeline returned a catalog of over 200,000 candidate 
variable objects, for which only very modest assumptions 
have been made about the variability type. This sample was 
further searched for objects which somewhere in the light 
curve showed no significant variations in a window contain- 
ing about half of all photometric points, but had one or two 
brightening episodes (to allow for binaries, which often have 



two peaks). Wozniak et al. (2001) describes the details of 
the selection process. Briefly, a brightening episode is de- 
fined as 3 or 4 consecutive points deviating respectively by 
4 or 3 (T upwards with respect to the baseline flux (deter- 
mined from the quiet period). At this point there arc no 
assumptions made about the shape of the light curve near 
the event. In all, 4424 light curves passed the above criteria. 
Further cuts requiring the light curves to be satisfactorily fit 
by the standard microlensing model would strongly discrim- 
inate against the non-standard ones such as parallax events 
(see Section 1). To avoid this problem, the most efficient way 
to recover as many of these events as possible, including the 
interesting non-standard ones, is still a visual search. 512 
events were found in the course of visual inspection of all 
4424 light curves in the wealdy filtered sample, the largest 
set of microlensing light curves published so fax. Allowing 
for slightly larger number of searched fields, this is roughly 
a factor of 2 more than discovered by the standard photo- 
metric pipeline from essentially the same data (Udalski et al. 
2000) . In contrast to the standard pipeline, the error distri- 
bution from the difference image analysis follows a Gaussian 
distribution and the photometric error is reduced by a fac- 
tor of 2-3. These two properties further increase the chance 
of seeing departures from the standard microlensing model. 
The error bars were re-calibrated so as to enforce the P^r 
degree of freedom to be unity for the best-fitting model (see 
§4). 



3 SELECTION PROCEDURE 

We first fit each microlensing candidate with the standard 
single microlens model which is sufficient to describe most 

events. In this model, the (point) source, the lens and the 
observer arc all assumed to move with constant spatial ve- 
locities. The standard light curve, A{t) is given by (e.g. 
Paczyhski 1986): 



A{t) 
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where wo is the impact parameter (in units of the Einstein 
radius) and 



r{t) 



t-to 
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with to being the time of the closest approach (maximum 
magnification), te the Einstein radius projected onto the 
observer plane, v the lens transverse velocity relative to the 
observer-source line of sight, also projected onto the observer 
plane, and the Einstein radius crossing time ^ . The Ein- 
stein radius projected onto the observer plane is given by 



te = 



c2(l-x) ' 



(3) 



where M is the lens mass, the distance to the source and 
X = Di/Ds is the ratio of the distance to the lens and the 
distance to the source. Eqs. (1-3) shows the well-known lens 



t In this paper, we differentiate tE from the duration of a mi- 
crolensing event, which is defined at the end of this section. 
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degeneracy, i.e., from a measured t-E, one can not infer v, M 
and X uniquely even if the source distance is known. 

The flux difference obtained from difference image anal- 
ysis can be written as 

fit) = /l [Ait) - 1] + A/, (4) 

where /l is the baseline flux of the lensed star, and A/ is 
the difference between the total baseline flux [i.e. the flux 
of the unlensed source and the blended star(s), if present] 
and the flux of the reference image. Note that in general 
A/ does not have to be zero (or even positive). Therefore, 
to fit the /-band data with the standard model, we need 
five parameters, namely, /l. A/, uo,to, and ts- Best-fitting 
parameters (and their errors) are found by minimizing the 
usual using the MINUIT program in the CERN library§. 

We then proceeded to fit these light curves with a model 
which accounts for the parallax effect. To do this we need 
to describe the lens trajectory in the ecliptic plane. This 
requires two further parameters, namely the projected Ein- 
stein radius onto the observer plane, te, and an angle tp 
in the ecliptic plane, which is defined as the angle between 
the heliocentric ecliptic x-axis and the normal to the tra- 
jectory (This geometry is illustrated in Fig. 5 of Soszynski 
et al. 2001). Once these two parameters are specified, the 
resulting lens trajectory in the ecliptic plane completely de- 
termines the separation between the lens and the observer 
(i.e., the quantity which is analogous to the standard model's 
Wo parameter from eq. 1). This allows the light curve to be 
calculated; the complete prescription is given in Soszynski et 
al. (2001), to which we refer the reader for further technical 
details (see also Alcock et al. 1995; Dominik 1998). For the 
parameters /l, A/, uo,to, and ts, we take the fit parameters 
from the standard fit as the initial guesses, while tb and 
are arbitrarily chosen (see below for more details). The best- 
fitting model is again found by minimizing the x^- While the 
standard fit to an observed microlensing light curve is almost 
always unique, this is not necessarily the case for the paral- 
lax events. To avoid missing the best-fitting parallax fit in 
the multi-dimensional parameter space, we ran the parallax 
fitting program with 24 combinations of te and ip as initial 
guesses; the model with the lowest is selected as the best 
parallax fit. 

We then compared the x^ values for the standard and 
parallax fits to determine which events were better described 
by the parallax model. Since the parallax fit utilized two 
additional parameters, we had to establish whether a given 
improvement in x^ was simply due to the increase in the 
number of free parameters, or whether this improvement was 
actually due to a deficiency in the standard model. This was 
done by performing the standard F-test for the significance 
of parameters (e.g. Lupton 1993) on each event. Briefly, we 
first calculate the variable 

[xi/Xp-l] X (iV-7)/2, (5) 

where xi and xp arc the for the best standard and par- 
allax fits, respectively, — 7 is the number of degrees of 
freedom for the parallax fit, and the factor of 2 is the differ- 
ence in the degrees of freedom between the parallax fit and 

§ http://wwwinfo.cern.ch/asd/cernlib/ 



the standard fit. The variable / follows an F2,n-7 distribu- 
tion. The statistical significance can then be evaluated using 
a one-sided test, i.e., the upper-tail of the F2,n-7 distribu- 
tion. A large value of I, i.e., a small one-sided probability, 
Pf = Fi> I), means that there is a significant improvement 
using the parallax model and vice versa. We selected events 
with probability less than pp < 0.001 as potential parallax 
candidates. In total, 109 microlensing candidates passed this 
selection criterion. We visually examined many of these can- 
didates, and found that the database was still contaminated 
by events with low signal-to-noise ratios. We needed to nar- 
row down the parallax candidates further and to do this we 
employed two additional criteria: the event 'duration' and 
the peak signal-to-noisc ratio. 

The duration constraint is imposed because the paral- 
lax effect should, in general, be more prominent for events 
that have long duration, since the Earth is able to move 
a substantial distance around the Sun during these events. 
However, the Einstein radius crossing time, ts, can some- 
times be misleading. For example, bright events with large 
peaJi fiuxes and small photometric error bars exhibit notice- 
able magnification for a period longer than ~ 2tE because 
the microlensing variability can still be detected even out- 
side the Einstein ring. To avoid this problem we clEissified the 
duration of an event as being the length of time which the 
flux of the standard fit was 3cr above its baseline value. This 
corresponds to the length of time which we can practicably 
utilize in our analysis of the Icnscd part of the light curve. 
For every event we recorded the duration along with the 
previously mentioned probability, and events with duration 
greater than 100 days were considered as potential parallax 
affected candidates. In total, 18 microlensing events passed 
this additional test. 

However, more than two thirds of these events still turn 
out to have noisy data and are unsuitable for detecting par- 
allax signatures. A typical example of this is given in Fig. 
1(a), which clearly illustrates the fact that more accurate 
data are needed in order to observe the often slight devia- 
tions between the standard and parallax fits. Another com- 
mon problem which occurs with the noisy events is shown in 
Fig. 1(b). In this case, the parallax fit substantially improves 
the x^ because the parallEix fit has multiple peaks that better 
reproduce the fluctuations in the data. However, one notices 
that the secondary peak falls in a gap between two observing 
seasons. Many other parallax flts for the noisy events show 
similar multiple-peak structures. While such multi-peak par- 
allax events have been predicted (Gould et al. 1994), the 
large error bars make the identification un-convincing. This 
is further supported by their often unusual parameter val- 
ues; for example, the majority of these events have values of 
te which are significantly different from the value of a few 
AU which is expected for a typical microlensing event. It is 
therefore clear that the identification of parallax signatures 
requires accurate data in order to observe the slight devi- 
ations and so the light curves for these noisy events were 
of little use and we do not regard them as parallax events. 
To eliminate these events from our selection procedure we 
needed to quantify how noisy the data were for an event. 
This was done by analysing the signal-to-noise ratio. To as- 
sess this signal-to-noise ratio, we use a quantity, (S/N)peak, 
defined as the difference between the peak flux and the base- 
line flux divided by the average error estimated using the 
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data points outside the 'duration' of the hght curve. We 
find empiricaUy that (S/N)peak > 30 provides an excellent 
separator. This criterion, while satisfactory, is somewhat ar- 
bitrary; we return to this issue briefly in the discussion. 

As well as identifying events which displayed paral- 
la^x signatures, we also wished to find events which unex- 
pectedly exhibited no deviations from the standard light 
curve. Of particular interest are events which have unusu- 
ally large time-scales but no signs of asymmetry, since they 
provide lower-limits on te (see §4.3). To isolate those we 
used the event duration which weis mentioned earlier. We 
also used two further criteria: firstly, a measure of the im- 
provement in between the standard and the parallax fits 
(i.e. a high value of pp, which corresponds to no signifi- 
cant improvement in x^) and secondly, the signal-to-noisc 
ratio, (S/N)peak. Events which had duration > 100 days, 
(S/N)peak > 30, and probability pp > 0.05 were considered 
as long-duration events that showed no parallax signatures. 



4 RESULTS 

Out of the 512 events, there were 8 events whore the incor- 
poration of the parallax efi^ect substantially improved the 
goodness of the fit (see above) and passed our test of event 
duration and the peak signal-to-noise ratio. These events are 
given in Table 1. One light curve was found to show clear 
signs of parallax affected behaviour, in addition to a num- 
ber of others which could be classified as marginal cases. 
In the following, we shall first discuss these parallax can- 
didates and in §4.3 we discuss the physical limits that one 
can derive from the long-duration events that exhibited no 
parallax signatures. 

For all the microlcnsing events presented here, we re- 
normalize the per degree of freedom to be unity for the 
best-fitting model by multiplying the quoted observational 
errors by a constant factor (usually between 0.8 to 1.3). 
This is necessary because we found that the quoted error 
bars in observations are often too large. The events in §4.1 
and §4.2 were re-normalized using the best-fitting parallax 
model, while the events in §4.3 were re-normalized using the 
best-fitting standard model (which has essentially the same 
X^ as the best- fitting parallax model) . This re-normalization 
affects primarily the quoted error bars in the model param- 
eters. It also weakly affects the selection procedure through 
both the peak flux ratio, (S/N)peak, and duration, although 
this perturbation is only very slight since I in eq. (5) remains 
unchanged. 

4.1 Parallax microlensing candidate 

From the catalogue of 512 microlensing events, the most 
convincing parallax event we found was sc33_4505, which is 
shown in Fig. 2. Despite the lack of data for two impor- 
tant regions (namely t = JD — 2450000 < 550 days and 
760 < t < 870 days), the systematic deviations from the 
standard curve can clearly be seen. During the upward slope 
of the light curve the standard fit undoubtedly has a gradi- 
ent which is too shallow and yet for the downward slope this 
fit's gradient is clearly too steep. This prominent asymmetry 
provides strong evidence that the light curve may be exhibit- 
ing parallax behaviour and this is further supported when 



the parallax model is used to fit the data. This model pro- 
vides a much more accurate fit, with the being drastically 
reduced from 398.972 (or x^ = 2.23 per degree of freedom) 
to 177 (x^ = 1 per degree of freedom). When observing re- 
sumed after t = 870 days it can be seen that the parallax 
model fitted this region more accurately, even though the 
lensing event was very nearly over and the magnification 
was only a fraction of the peak value. 

The best fit parameters from both the standard and 
parallax fit are given in Table 2. In particular, for the best- 
fitting parallax model, we find that 

tE = 194 ±30 days, = 6.37 ± 0.53 AU. (6) 

An expression for the lens mass as a function of the relative 
lens-source distance can be shown to be (see Soszynski et al. 
2001), 

«=#(^-i)=— (S)'(oS;).- 

where TTroi is the relative lens-source parallax in milli- 
arcseconds. For this event, the parallax parameters give a 
projected velocity of 

v = ^= 57t\t km s-\ (8) 

IB 

The lens mass depends on the relative lens-source parallax 
(eq. 7). If the source is about 7kpc away, and the lens lies 

in the disk half-way between the observer and the source 
{x = 1/2), then TTrd ~ 0.143 mas, which gives a lens mass 
of about O.71Af0; as a comparison, for a bulge self- lensing 
event with Ds w 8kpc and Di w 6kpc, then TTrei 0.042 
mas, which would give a lens mass of about O.21M0. How- 
ever, this latter scenario is less likely since the projected 
velocity of the lens is relatively low. It appears that similar 
to the previous parallax microlensing events, this event is 
caused by a slow-moving and likely low-mass lens. 

4.2 Marginal parallax microlensing candidates 

Apart from the above convincing candidate, sc33_4505, a 
number of events wore found which showed indications of 
parallax signatures. However, the case for clEissifying these 
events as having parallax affected light curves is not wholly 
convincing, mainly due to the insufficiency of the data. Figs. 
3-9 show these 'marginal' microlensing events that we have 
found, and Table 3 presents the best fit parameters for both 
the standard and parallax models. 

The problem of insufficient sampling is clear from sev- 
eral marginal cases. For example, for the sc41_3299 event 
(Fig. 8), there are very few points in the rising branch 
of the light curve. As a result, although the parallax fit 
seems to provide a somewhat better fit in the region be- 
tween t = 880 — 900 days, it is not clear whether this is 
a genuine parallax event or just due to the largo errors in 
the data. The situation is similar for events sc20_5748 (Fig. 
4), sc27_3078 (Fig. 6), and sc43_836 (Fig. 9). For the event 
sc6_2563 (Fig. 3) , the standard fit consistently over- predicts 
the fluxes between t = 1040 to 1140 days, and the parallax 
model provides a much better fit for those parts. Unfortu- 
nately, there was no data for the time period {t ~ 1180 days) 
when the standard model and the parallax model are pre- 
dicted to show substantial deviations. The fact that we can 
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pick out marginal cases like sc41_3299 implies that our se- 
lection method is sensitive, and it is unlikely that we could 
have missed many true candidate parallax events. 

The situations for the events sc20_5875 (Fig. 5) and 
sc35_2526 (Fig. 7) arc slightly different. For the sc35_2526 
event, the light curve is highly asymmetric, and the stan- 
dard model is clearly a bad fit. The parallax model, on the 
other hand, provides a much better fit; the improvement in 
is from 361.5 to 170. However, the fit is not perfect, par- 
ticularly at the peak and for the few data points between 
700 < t < 750 days. For this reason, we still cautiously re- 
gard this event as a "marginal" parallax candidate event. 
The situation for sc20_5875 is similar, bxit less dramatic. It 
is possible that these two events arc produced by binaries, 
where the asymmetry is provided, e.g., by the shear per- 
turbation by the secondary on the primary lens (e.g. Mao 
& Di Stefano 1995). Notice that, for the event sc35_2526, 
the parameter is unusually long (see Table 3); however, 
this is just an artifact due to the well-known degeneracy 
between the impact parameter, uq, and for microlensing 
light curves in blended light curves (WozniaJj & Paczyhski 
1997), where only the combination of Mote can be inferred 
in some cases. 

4.3 Long microlensing events that show no 
parallax signatures 

Of almost equal interest to parallax events are the events 
which have large time-scales and yet exhibit no signs of par- 
allax induced asymmetry. A total of 20 events were identi- 
fied using the selection criteria which wore given in §3, i.e. 
events which had duration > 100 days, (S/N)peak > 30, and 
probability pr > 0.05. The parallax signatures become more 
pronounced for small values of fs because the Earth's mo- 
tion around the Sun then becomes a larger fraction of the 
Einstein radius, and hence presents a larger perturbation 
to the light curve. The lack of parallEix signatures therefore 
provides a lower- limit on tb. For each such event, we fix 
in a range of values while letting all the other six parameters 
vary and find the minimum x'^ value for each fixed rs- The 
2fT confidence constraint on fe is then given by the te value 
corresponding to the point at which its becomes larger 
than the best fit by 4.00 (e.g. Lupton 1993). From eq. (7), 
one sees that this lower limit on te can then be translated 
into a lower limit on M/niei- The events which produced 
the 4 largest constraints are given in Table 4, and their light 
curves are shown in Fig. 10. 

In all cases, the lower limits on the mass are rather 
weak, the best case being for event sc26_2218, whore the 
lower limit reaches about 0.2Mq for TTrei = 0.1 mas. These 
weak constraints are similar to those found by the EROS col- 
laboration for the microlensing event ER0S2-GSA1 (Derue 
et al. 1999), from which they obtained values of te > 1-33 
AU and IO-^Mq. 



5 SUMMARY AND DISCUSSION 

We have conducted an extensive search for parallax events 
in the microlensing database of 512 microlensing events ob- 
tained by the difference image analysis of the OGLE II data. 
Our selection procedure involves a direct comparison of the 



standard microlensing model and the parallax model, aug- 
mented by the requirements of long duration and high peak 
signal-to-noise ratio. While empirically we found these cri- 
teria to be very effective, there is no guarantee that these 
selection procedures are optimal. This issue is best explored 
using the recovery of artificial parallax events in Monte Carlo 
simulations (see below). 

Using our selection criteria, we found one convincing 
new parallax microlensing event, sc33_4505. A number of 
other more marginal cases were also recovered. However, 
while these marginal events are better fitted by a model 
incorporating the parallax effect, their parallax nature can- 
not be established beyond any doubt due to the poor sam- 
pling in the light curve. This may be particularly severe for 
the events which last for about one year, due to the gap 
of approximately three months between observing seasons. 
In such cases either the rising or the declining branch of 
the light curve will typically be missing, and so the paral- 
lax induced asymmetry is particularly difficult to identify. 
We have also found some very long events that showed no 
parallax signatures, and these events provide a lower-limit 
on the Einstein radius projected onto the observer plane. 
Consequently, a lower-limit on the lens mass M /vTrei can be 
derived (cf. eq. 7). However, in most cases, the lower limits 
on M/TTrei are rather weak. 

The convincing parallax event (sc33_4505) is similar to 
all of the other 4 known parallax events (see introduction) 
in that it is caused by a slow-moving lens and, quite likely, 
a low mass lens. The low projected velocity favors the inter- 
pretation that this microlensing geometry is a disk source 
lensed by a disk lens. For such events, the observer, the lens 
and the source rotate about the Galactic center with roughly 
the same velocity, and the relative motion is only due to the 
small, ~ lOkms"^, random velocities (see, e.g. Derue et 
al. 1999). On the other hand, the chance for a bulge source 
(with its much larger random velocity, ~ 100 km s~^) to have 
such a low projected velocity relative to the lens (whether in 
the disk or bulge) is small. So it appears that all known par- 
allax events axe caused by disk-disk lensing. If this is true, 
then the radial velocities of the lensed sources are expected 
to be small. It will therefore be very interesting to check this 
by obtaining the radial velocities of the parallax events spec- 
troscopically. Since their projected transverse velocities are 
known, one obtains a more complete kinematic picture of 
these unusual microlensing events that can be used to probe 
the dynamical model of the Milky Way. As a by product, 
one also obtains the metallicity and age of these stars. 

While the number of convincing parallax events toward 
the Galactic bulge seems to be very low (1 out of 512), the 
number of marginal cases makes the true fraction somewhat 
uncertain. A related question is whether existing microlens- 
ing catalogs are biased against parallax events. This is an 
important issue because parallax events preferentially have 
long durations, which may be detectable only when one has 
monitored the stars for many years. In fact, after the com- 
pletion of this work, eight additional events were recovered 
by cross identification of the 214 microlensing candidates 
from the standard OGLE database with the difference im- 
age analysis variability database of Wozniak (2000). Some of 
these were missed in the first search because they have such 
long durations that they have not yet reached a constant 
baseline. We plan to analyze these additional unique events 
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for parallax signatures. On the theoretical side, Buchalter & 
Kamionkowski (1997) have estimated the expected fraction 
of parallax events which would be identified, using a regu- 
lar sampling of every few minutes to ~ 1 day. Our search 
strongly suggests that the irregular sampling and gaps be- 
tween different observing seasons hamper the recovery of 
parallax events. We plan to perform a simulation with real- 
istic sampling and photometric errors similar to those found 
in observations. A comparison between the observed and 
predicted rates may enable us to provide constraints on the 
lens and source kinematics. Such simulations will also be 
helpful for devising the optimal search strategy for selecting 
parallax events from the experiments. The results of these 
Monte Carlo simulations will be reported in a future publi- 
cation. 
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sc33_4505 
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Figure 1. Two events which highlight the problems associated 
with noisy data. The solid and dotted lines represent the stan- 
dard fit and the parallax fit, respectively, (a) The event sc39_5310, 
which illustrates how slight deviations between the two fits can 
easily be hidden by noisy data. The error bars have been omitted 
for clarity, but a typical error bar is shown in the top left hand 
corner, (b) The event sc30_3395, which demonstrates how a mi- 
crolensing light-curve can be seemingly-erroneously fitted with a 
parallax model exhibiting many peaks. This leads to a significant 
improvement in x^- Such multi-peak parallax models often occur 
in events with noisy data, which may signal additional variability 
other than microlensing. 



Figure 2. Difference flux (/(t), as defined in eq. 4) as a func- 
tion of time for the best candidate parallax microlensing event, 
sc33_4505. The solid line shows the best standard fit while the 
dotted line shows the best fit that accounts for the parallax ef- 
fect. The insets show the light curves for the 'non-lensed' seasons. 
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Figure 3. The marginal candidate parallax microlensing event, 
sc6_2563. The notations are the same as in Fig. 2. 
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Figure 4. The marginal candidate parallax microlensing event, 
sc20_5748. The notations are the same as in Fig. 2. 



Figure 6. The marginal candidate parallax microlensing event, 
sc27_3078. The notations are the same as in Fig. 2. 
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Figure 5. The marginal candidate parallax microlensing event, 
sc20_5875. The notations are the same as in Fig. 2. This event is 
not perfectly fit by the parallax model. This event may, instead, 
be caused by a (weak) binary lens (see also Fig. 7). 



Figure 7. The marginal candidate parallax microlensing event, 
sc35_2526. The notations are the same as in Fig. 2. This event is 
not perfectly fit by the parallax model. This event may, instead, 
be caused by a (weak) binary lens (see also Fig. 5). 
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sc41_3299 




Figure 8. The marginal candidate parallax microlensing event, 
sc41_3299. The notations are the same as in Fig. 2. 

sc43_836 




Figure 10. The 4 long-duration lensing events which produced 
the most stringent lower-limits on the lens mass M/n^^i (see eq. 
7 and Table 4). For these events, the incorporation of parallax 
effects produces virtually no improvement in the x^- The dotted 
line shows the theoretical baseline flux and the horizontal axis 
covers a length of time corresponding to the event duration, as 
specified in §3. 
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Figure 9. The marginal candidate parallax microlensing event, 
sc43_836. The notations are the same as in Fig. 2. 
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Table 1. The 8 events which were systematically selected as po- 
tential parallax affeeted light-curves, and their subsequent classi- 
fication. The probability pp and (S/N)p(,ak a-re defined in §3. The 
convincing candidate (sc33_4505) and other marginal candidates 
are shown in Figs. 2-9. 



Event RA DEC Duration (day) tg (day) (S/^)pcak Classification 



sc6_2563 18:08:02.64 


-31:49:05.2 245 66.6 ± 1.3 1.22 X 10" 


-8 


228.39 


marginal 


sc20_5748 17:59:39.05 


-28:27:16.6 150 78.2 ± 2.5 2,58 X 10" 




38.99 


marginal 


Bc20_5875 17:59:08.99 


-28:24:54.7 146 27.22 ± 0.21 0.00 




119.04 


marginal /binary? 


sc27_3078 17:48:15.92 


-34:50:09.0 162 124 ± 10 2.04 X 10" 


-7 


142.82 


marginal 


sc33_4505 18:05:46.71 


-28:25:32.1 663 194 ± 30 5.78 X 10" 


-32 


78.79 


convincing 


sc35_2526 18:04:22.42 


-27:57:52.2 162 6600 ± 180 1.44 X 10" 


-28 


47.02 


marginal /binary? 


8C41_3299 17:52:19.08 


-32:48:20.8 119 98 ± 18 1.28 X 10" 


-5 


35.17 


marginal 


sc43_836 17:34:52.64 


-27:23:14.4 207 45.1 ± 1.6 4.31 X 10" 


-5 


31.70 


marginal 


Table 2. The best standard model (first row) and the best par- 








allax model (second row) for the convincing candidate sc33_4505. 








Model tQ 


*E (day) ug /l A/ 






fE (AU) 


S 647.22 ± 0.36 


166 ± 11 -0.412 ± 0.036 730 ± 91 -172.2 ±1.8 






— 399.0 


P 654.4 ±4.1 


194 ± 30 -0.207 ± 0.065 440 ± 150 -174.3 ±2.1 


3.136 ± 0.011 


6.37 ±0.53 177 



Table 3. The best standard model (indicated by 'S') and the 
best parallax model (indicated by 'P') for marginal parallax can- 
didates. The error bars have been rescaled so as to enforce the 
per degree of freedom to be unity for the best-fitting parallax 
model. 



Event 


to 


*E (day) 


UQ 


/l 


A/ 








fE (AU) 




sc6_2563,S 


1251.971 ± 0.068 


89-2 ± 2.8 


0.0758 ± 0.0032 


135.9 ± 5.7 


-49.43 ± 0.52 










253.7 


P 


1263.85 ± 0.62 


66.6 ± 1.3 


0.0272 ± 0.0085 


192.8 ± 5.5 


-48.84 ± 0.49 


0.262 ± 


0. 


,012 


4.27 ± 0.16 


214 


sc20_5748,S 


1277-03 ± 0.22 


357 ±120 


-0.0148 ± 0.0056 


4.1 ± 1.5 


-26.24 ± 0.67 










243,2 


p 


1315.12 ± 0.96 


78.2 ± 2.5 


-0.325 ± 0.041 


16.67 ± 0.79 


-25.50 ± 0.41 


1.234 ± 


0. 


.057 


2,43 ± 0,11 


221 


sc20_5875,S 


948.74 ± 0.10 


31.3 ± 1.7 


0.461 ± 0.042 


1950 ± 250 


-181.9 ±1.7 










666,7 


P 


933.88 ± 0.28 


27.22 ± 0.21 


2.073 ± 0.011 


320 ± 25 


-182.0 ±1.6 


0.235 ± 


0. 


Oil 


0.3515 ± 0.0021 


189 


sc27_3n78,S 


1243-481 ± 0.077 


312 ± 77 


0.0070 ± 0.0018 


4.8 ± 1.2 


-12,08 ± 0.52 










245,2 


p 


1238-75 ± 0.34 


124 ± 10 


-0.6266 ± 0.0068 


18.2 ± 1.3 


-11.60 ± 0.36 


2.971 ± 


0. 


,026 


1.68 ± 0.019 


212 


sc35_2526,S 


696.33 ± 0.20 


53000 ± 11000 


0.000107 ± 0.000021 


0.0291 ± 0.0057 


-7.39 ± 0.59 










361,5 


p 


695.46 ± 0.23 


6600 ± 180 


0.01512 ± 0.00018 


0.0920 ± 0.0066 


-4.65 ± 0.57 


2,9967 ± 


0. 


.0087 


63,02 ± 0,73 


170 


sc41_3299,S 


927.65 ± 0.29 


63.2 ±8.3 


0.186 ± 0.034 


40.8 ± 8.5 


-4.63 ± 0.40 










197,1 


P 


913.1 ± 7.7 


98 ± 18 


0.19 ± 0.16 


26.2 ± 6.8 


-4.87 ± 0.47 


6.222 ± 


0. 


.096 


2.33 ± 0.68 


173 


sc43_836,S 


1054.69 ± 0.41 


71.5 ± 7.2 


0.437 ± 0.067 


126 ± 28 


-17.91 ± 0.41 










262.0 


,P 


1040.0 ± 1.6 


45.1 ± 1.6 


-0.104 ± 0.035 


176 ± 13 


-17.38 ± 0.37 


3.1744 ± 


0. 


,0054 


1.511 ± 0.091 


241 



Table 4. The 4 non-lensed events which produced the best 2a 
confidence lower limit on the projected Einstein radius in the 
observer plane (fg) and on M/ir^^i (cf. eq. 7). The 'duration' 
of the event is defined in §3, and the probability pp is defined 
below eq. (5). The relative lens-source parallax TTrei is in units 
of milli-arcseconds; for microlensing events toward the Galactic 
bulge, TT^ei ~ 0.1 mas. 



Event 


RA 


DEC 


Duration (day) 


tE (day) 


PF 


'■E 


(AU) 


"/"rel (M©) 


sc22_1992 


17:56:35.30 


-30:56:33.4 


151 


45.5 ± 3.7 





.132 





.94 


0. 


11 


sc26_2218 


17:47:23.29 


-34:59:52.4 


195 


45.10 ± 0.59 





.104 


3 


.74 


1. 


71 


sc34_3906 


17:58:37.11 


-29:06:29.9 


331 


58.7 ± 6.2 





.204 


1 


.44 


0. 


25 


sc45_810 


18:03:30.06 


-30:09:55.6 


115 


98 ± 54 





.151 


1 


.73 


0. 


37 
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